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Normal modes and their propagation constants have been found for a 
two-dimensional lens-like medium in which the transverse refractive index 
varies essentially with quadratic law but has perturbing terms of higher 
order. In such a realistic guiding medium with aberrations, those modes 
are used to find the field configuration of a Gaussian beam of half width W 
entering off-axis. 

Close to the input the beam oscillates periodically with amplitude x f 
as if the medium were aberration-free, but slowly the beam cross-section 
changes shape, breaking up in several maxima, and increases size, reaching 
a maximum approximately equal to 2(W -f- x { ). Afterwards the beam 
slowly shrinks back to the starting field configuration and the process 
repeats again. 

These results are applicable to a sequence of lenses ivith aberrations 
and become important when the lenses are closely spaced. If redirectors 
are to be used to compensate for lens misalignments, the corrections must 
be made before a large break-up of the beam occurs. 

I. INTRODUCTION 

Transmission through a leus-liko medium consisting of a dielectric 
rod in which the dielectric constant decreases radially with quadratic 
law has been studied because, first, it is closely related to light transmis- 
sion in a periodic sequence of gaseous lenses, 1,2 ' 3 second, it helps to 
understand the filamentary nature of the oscillations in ruby lasers, 4 
and third, if the dielectric constant is complex, it describes the gain 
medium in a gaseous laser. 5 

Whenever the radial dependence of the dielectric constant is quadratic, 
any paraxial beam propagates undulating periodically about the axis, 
and the field reproduces itself after each period. This is true only within 
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some approximations.* Part of the object of this paper is to find the 
distance over which the approximations hold and how it depends upon 
the beam shape and displacement from the axis. The mathematical 
description of Gaussian beams in such an aberration-free lens-like 
quadratic medium has been achieved via two techniques. The first 
consists of expanding the input beam in terms of on-axis normal modes 1 
and summing them up everywhere. The second technique consists of 
solving differential equations that elegantly relate the location of the 
beam axis, the spot size, and the curvature of the phase front at any 
place along the beam. 7 Nevertheless, we show in this paper that if the 
lens-like medium has aberrations, no matter how small, the electro- 
magnetic field of a paraxial beam changes shape radically as it travels 
along the guide. The distorted beam cannot be described any more by 
the beam axis location, the spot size, and the phase front curvature, 
and consequently, the off-axis beam must be calculated following the 
first technique. 

The on-axis normal modes and their propagation constants are found 
in Section II using a first-order perturbation technique; the off-axis 
beam is calculated in Section III; the limits of validity of the previous 
results are given in Section IV and conclusions are reached in Section V. 

The modes in lens-like media with arbitrarily large perturbations 
have been studied by S. E. Miller 8 and J. P. Gordon. 9 However, they 
have not given a quantitative description of an off-axis beam. 

Beam deformations similar to those analyzed in this paper have 
been obtained by D. Marcuse 10 who calculated, via a computer, the 
transmission of an off-axis Gaussian beam through a sequence of 
curved and distorted lenses. 

II. MODES IN A LENS-LIKE MEDIUM WITH SMALL ABERRATION 

Consider a lens-like medium in which the refractive index 



n = n„ 1 — 



,[i-(f)'~ L «.(?)' 



(1) 



* A. H. Carter has found 6 that if the dielectric constant of the lens-like medium 
varies radially, not with quadratic law, but with the square of the hyperbolic secant 
law, the periodicity is not restricted to paraxial beams and since there are no approxi- 
mations involved, the results hold for any length. Unfortunately, the functions de- 
scribing the modes are hypergeometric and they cannot be easily handled as the 
parabolic cylinder functions of the quadratic medium. Because of this and also be- 
cause our main interest is the study of transmission through non-ideal media anyhow, 
we will keep on basing our calculations on the quadratic medium and we will even 
call it ideal. 
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varies only in the transverse direction x; n is the refractive index 
at x = 0, and L and a a are constants whose physical significance will 
be shown later. The quadratic law of the dielectric constant is slightly 
perturbed by the summation in the integers a = 1, 2, 3 • • • . 

We study here modes propagating along z in a medium with variation 
only in the x direction; nevertheless, the extension to the case in which 
the dielectric constant varies also along the other transverse direction y, 
can be easily achieved as shown in Ref. 1. 

Assuming that the only component of electric field E is parallel 
to the y axis and independent of y the wave equation is 

dx 2 + dz 2 + U 



1 ~ T ~ Z^aA-r 



L J ^ " Q VL 



E = 0, (2) 



where X is the plane wave wavelength in a medium of refractive 
index n . 
Let us call 

f = | (3) 

a new transverse variable normalized to the beam half-width 

in the unperturbed medium, and also 

E = E p exp(-iy^). (5) 

E„ , the transverse field distribution of the pth mode, is a function 
of £ exclusively and 

y.-*yl 1 -( p + i+ Za J .®"-')$ (6) 



is the propagation constant of the pth mode along the z axis. The 
summation ^« a a f a (\/L) a/2 ~ l , is the perturbation due to the aberra- 
tions and the unknown /„ will be determined later. 

Now we substitute (3), (4), and (5) in (2) and obtain an equation 
in | exclusively for the transverse field E p , 



f + { f+ !-(!)- +? ..Q 



(I)'] 



/. -U #p = 0. (7) 



For a a = 0, the solution 11 is the parabolic cylinder function Z) p (£) = 
ex P ( — £ 2 /4)# p (£) which is a product of the Gaussian function exp ( — £ z /4) 
and the Hermite polynomial of integer and positive order p. 
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For a a ^ it is shown in Appendix A that if 

z^p) b/2, «i (8) 

the differential equation (7) can be solved by using first-order stationary 
perturbation theory. 12 Then, the transverse field distribution of the 
pth mode in the perturbed medium expressed in terms of the normal 
modes of the unperturbed medium is 

E P = D,(©+ Zc„D n (S); (9) 

ny*v 

the sum extends from to «> excluding n = p and the values of c n 
are given in (54). We also find in (55) that 



/. = - 



2~* a/2 a\ p! 2 / \ if a is even » 

"■°(m!) a (p-m)!(|- m)\ (">) 



if a is odd. 



Further on, specific values of /„ will be needed, so some of them are 
given in Table I. From the table and (6) we conclude that if the medium 
has only antisymmetric perturbation (a odd), the propagation con- 
stants of its modes are identical to those of the modes in the unperturbed 
medium. Physically, this means that to a first order, the change in 
phase velocity introduced by the perturbation on one side of the guide 
is canceled by the other. 

III. BEAM IN THE PERTURBED MEDIUM 

We want to calculate the field F{z) everywhere in the perturbed 
medium such that at the origin it coincides with a prescribed function 



Table I- 


- Some Specific Values of /„ 


a 


fa 


1 





2 


i(l + 2p) 


:>, 





4 


A(l + 2p + 2p 2 ) 


6 





6 


4f(l+ fp + 2p 2 + *P 3 ) 
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F(0). In terms of the normal modes E p , the field is 

F{z) = Z-lp^expC-tV). (11) 

p-0 

and at the origin 

F(0) = ± AJE V . (12) 

p-0 

We simplify next A P E P and the propagation constant y p . 

Since the approximately calculated normal modes of the medium 
with aberration are orthogonal to the first order of the perturbations 
a a (\/L) a/2 ~ l , we calculate, with the help of (9), each term of the 
summation to be 

A P E P = B P DM + 0[a a (\/L) a/2 - 1 }, (13) 

where 

f now,® # 

^ = J -^—~ (14) 

V — CO 

is the amplitude of the pth mode of an expansion of the input field 
F(0) in terms of modes of the unperturbed guide and the term of the 
order of a„(\/L) a/2 ~\ given only for completeness, is 







a "\I 



XV 



- D.® f F(0)D p (it)dt+D p (i;) f F(0)D n ®di: 






Substituting (13) in (11) one obtains 

m = Z BJ>JQ exp (-iyj) + £ o[a„(^) n/ 
p-o p=o L v/-" 



exp(-zV)- (15) 



We will find that for large z the first summation yields configurations 
which depart grossly from the field in the unperturbed medium, there- 
fore, according to (8), the second summation whose amplitude is of 
the order of a a (X/L)° /2 ~ l will be neglected. 

Furthermore, the propagation constant y P , (6), which is exact for 
the ideal quadratic medium and good to first order of ^ Q a„(\/L) a/2 ~ l 
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can be expanded in series 



X 



![ P+ i + ?../.sri£- •}• 



(16) 



The third and higher-order terms can be neglected only if their con- 
tribution to the phaseshift of the highest-order mode P of significant 
amplitude is small, that is, if 



tt\z 
4L 2 



1 . t- . ( X 



p + l+IX/.u 



a/2-1 



«1. 



According to (47) in the Appendix, once the inequality (8) is satisfied, 
the summation is negligible compared to P and consequently the in- 
equality above reduces to 

jf- 2 P 2 «\. (17) 

4L 

This expression establishes the range of validity z, of the description 
of a beam composed essentially of P modes at wavelength X, traveling 
in a lens-like medium characterized by L. Since the inequality is in- 
dependent of a a we conclude that the finite range of validity z is deter- 
mined not by the perturbation of the medium, but by the approximation 
involved in the expansion (16) of the propagation constants thus 
affecting the description of the beam even in the ideal quadratic medium. 
P Once P is known, the inequality (8) determines the magnitude of 
the perturbation a a for which the calculations are valid. Then, provided 
that inequalities (8) and (17) are satisfied the simplified version of the 
field (15) becomes 



F(z) = Z B P D V (Z) exp \-\yvi ~ I Z <*«/«(£) 



(18) 



in which 



2t 



1 -£&> + *>] (I") 



is the propagation constant of the pth mode in the unperturbed ideal 
quadratic medium. 

Therefore, the field F(z) is described by a summation of modes of the 
ideal quadratic law medium but Avith phase constants corresponding to 
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modes in the perturbed medium. Actually, if a a = 0, we derive from 
(18) the field in the ideal quadratic medium 

I® = Z BA<8 exp (~iy pi z). (20) 

In order to find out more about the field distribution along the 
perturbed guide, we will work out a typical example assuming that the 
only perturbation of the dielectric is of order 4. Dropping the subscript 
from a 4 and with the help of Table I we find that 

Ea4|P'=a/ ) | = |a|(l + 2 P + 2 P '). 

The field (18) then results 

F{z) = exp (i | ~ 

where the distance 



(•' i 7)) t BJ>& exp {-{ 7p , -~(P + p 8 )}} , (21) 



has a physical significance to be described later. 

I do not know how to add (21) in general but the summation can 
be performed for discrete and significant values of z. They are 

*„., - (m + 2")D (23) 

where n is an arbitrary integer and v is a positive integer. The summation 
can also be performed for z = (p — 2 - *).D but the results are quite 
similar to those found for z„, . For the particular values of z, given 
in (23), the exponential in (21) containing p 2 can be expressed as a 
sum of exponentials containing only p. As a matter of fact, 

exp [ivp 2 (ft + 2"")] 

= exp S —iirp 



M + 



(-")L 



\ 2-1 

\ £ G„ exp (-itf-'irgp) (24) 

J a-o 



and 

2'-l 



G„ = 2" 2 ^p [wr2"'s(s + 2q)\. (25) 

5-0 

The correctness of this expansion can be verified by substituting (25) 
in (24) and performing the summation first in q and then in s. 
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The field (21) at z = «„, results with the help of (24) 
Ffa,) = exp 



t|(M + 2-) 



where 



and 



EG. exp («/(<■„,.), ( 26 ) 

a-0 

r( 4L\fl_ _ _J SlI 

~ 2 V X/_2' (-")! 2-'J 



(27) 



U = (M + 2-1D + (L - ^1- - JL)L. (28) 

Ignoring the overall uninteresting phase 7t/2(m + 2""), equation (26) 
establishes that the field at z„, is made of 2" terms, each given by a 
factor G q exp i<p rQ times the field distribution /(£„„). As we saw before, 
this field (20) coincides with the field distribution that the input F(0) 
would have in the aberration free medium at a distance z = £„,„ given 
in (28). Therefore, in the perturbed medium, at z = z„, , the field 
F(z„„) due to an arbitrary input F(0) is described by the superposition 
of field distributions that the same input F(0) would produce in the 
unperturbed medium at 2" cross-sections located at distances f„»« from 
the origin, multiplied by certain phase shifts and amplitudes. 

Let us extend the example assuming the input F(0) to be an off-axis 
Gaussian of half-width W = VxL/t. In the aberration-free medium 
the beam trajectory is a sinusoid of period 2L, and at all cross-sections 
the field is Gaussian of half-width W (see Fig. 1(a)). Now we proceed 
to calculate the field F(z M „) in the perturbed medium at the specific 
abscissas z„, (23) assuming the same off-axis Gaussian input F(0). 

For v = 0, that is at abscissas 

2mo = G* + l)D (29) 

the field in the perturbed medium is derived from (26) and (25) to be 

FM = exp[t|(M+l)]/(f, uo )- (30) 

Except for the phase shift tt/2(m + 1) the field in the perturbed medium, 
F{z^) coincides with the Gaussian field in the unperturbed medium, 
Itt»oo) at the same abscissa 

*. = r„ uu - 0* + i)D. (3i) 

The distance D = 8L 2 /3aX is then the distance between successive 
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INPUT FIELD 

E - e"(w) 



t/e*) (<- 
< 



2W 



Sooo - z c 






r - D L 

— ?020 = ? + V— « 

C ^021- 4 4 * 

<" 5 022 --jj 4-*! 

r _ D 5i_ 

!"023 - 4 T" , 

e -.. 






(a) 

»7, 




-Id 

z 00 = D 

HORIZONTAL SCALE IS GROSSLY REDUCED 



Fig. 1 — (a) Off-axis beam in ideal medium, n = n B [l — (tx/l) 2 ]'. (b) Off-axis beam 
in perturbed medium, n = w„[l — (jtj/L) 2 — a(vx / L)*]i . 

D = 8L*/3oX and ?/> = VxT/x. 



Gaussian field distribution in the perturbed medium. We will call D 
a pseudo-period because the fields at abscissas z„ , in general, do not 
have the same position with respect to the axis as the input field. Only 
if 4L/3o\ is an integer those fields arc identical and the pseudo-period 
becomes a true period. 

Do avc know about the field close to z» , say within a few LI The 
perturbed medium differs only slightly from the quadratic one, there- 
fore close to z„„, the field in both media must be quite similar. The 
similarity of fields in both media in the neighborhood of abscissas 
(/x + 1)7) for n equal to —1 and is depicted in Figs. 1(a) and 1(b). 

For v = 1, that is at abscissas 



2 M1 = Gti + \)D 



(32) 
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half way between the previous ones (29), the field derived again from 
(26) and (25) is 

exp \i | 0* + 1) 
*•(*„) = L 



[i(r„io) - n(Ui)] 

(33) 



V2 

Ignoring the phases the field in the perturbed medium F(z„,) is made 
of the superposition of the fields I(?„ l0 ) and 7(f,„) in the unperturbed 
medium, reduced in amplitude by V 2 anfI found at abscissas derived 
from (28) to be 

U = (m + h)D + k (34) 



and 



= (M + h)D - \ (35) 



Again, since the medium is only slightly perturbed, the two Gaussians 
can be treated independently. Therefore, within a few L from z„, , 
the field in the perturbed medium consists of two Gaussian beams 
interleaving. In Fig. 1 (a), for ju = 0, the fields at abscissas f „io = (D+L)/2 
and fou = (D - L)/2 are located, and then used in Fig. 1(b) to con- 
struct the field in the perturbed medium around z 01 = D/2. 

For v = 2 we find that the field at z„ 2 = (ji + i)D, can be synthesized 
by the superposition of four fields that, except for the amplitude and 
phase, are found in the unperturbed medium at abscissas 

L 

_ L 

4 (36) 

3L 

S»i22 2»2 j » 

and 



5L 

f»<23 — Z(i2 A 

The field close to z„ 2 consists of four Gaussian beams weaving as shown 
in Fig. 1(b). 
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Within a pseudo-period the field in the perturbed medium differs 
substantially from the Gaussian input, and at a single cross-section 
it can exhibit several maximas if the peaks of the several Gaussians 
that depict the field at that cross-section are resolved. Furthermore, 
the cross-section of the beam becomes at most as wide as twice the 
half-width of the input beam plus twice the amplitude of the beam 
oscillation in the unperturbed medium. 

Even though the fields in the perturbed medium at distances close 
to D/2 and D/\, Fig. 1(b), are made of symmetrically interleaving 
beams it must not be concluded that the fields at D/8, D/16, etc. 
must be symmetric with respect to the x = plane. As soon as the 
interleaving beams overlap, the relative phases become important and 
the apparent symmetry breaks up. As a matter of fact, close to z = 0, 
the field (20) can be expressed in two ways. Either choosing n = — 1 
and v = in which case the field is that of a single beam or else picking 
n = and v — > <» f in which case the same field is made by the super- 
position of 2" interleaving beams with deceiving symmetry. 

If instead of Gaussian the input field has another shape such as that 
in Fig. 2(a), the field in the unperturbed medium reproduces the input 
at even multiples of L, while it repeats the input mirrored in the plane 




D/2- ~*j 

HORIZONTAL SCALE IS GROSSLY REDUCED 



Fig- 2 — (a) Off-axis beam in ideal medium, n = n„[l — (w/L) 1 ]!. (b) Off-axis beam 
in perturbed medium, n = n„[l - (xt/L) 2 - a(irx/Ly]K 
D = SL 2 /3aX and 4L/3a\ integer. 
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x = 0, at odd multiples of L. In between the beamwidth varies period- 
ically with period L and we represent it qualitatively in Fig. 2(a). 

The field in the perturbed medium at z <° , Zoo , «oi , and their neighbor- 
hood is depicted in Fig. 2(b) following similar steps to those taken 
for the previous example. 

If the perturbation of the medium were of sixth-order, instead of 
fourth, the pseudo-period, deduced from (18) and Table I, would be 
32/5 L 3 /a 6 \ 2 and as before, an off-axis beam would periodically deform 
itself exhibiting several maxima and the largest beam cross-section 
would be roughly twice the input beam half-width plus twice the 
input beam displacement from the axis. 

V. LIMITS OF APPLICABILITY 

Over what length z is the field in the unperturbed medium (20) 
valid? That length, calculated from the inequality (17), is 

2«^4- (37) 

We need the value of P, that is the highest-order mode of significant 
amplitude. 

Continuing with the example, the input 

F(G) = D (t - £.) 

is a Gaussian beam of half-width W, displaced x { (normalized £,) from 
the z axis. Then according to (14) the amplitude of the pth mode is 



B. 



/- exp [" 



ft - 6)' 

4 



l>:, </: /fcX . exp [ - | 



£ 



D'ft) 4 



(38) 



For 

p» 1 

the asymptotic value of B p results 13 in 

& 




(39) 

Obviously B v decreases rapidly for p such that the second parenthesis 
is smaller than one. Therefore, we select the highest-order mode of 
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significant amplitude to be 

p-f 

or substituting £,• by its equivalent in (3) 

p - i- (*» 

With this value and that for W given in (4) the range of validity (37) 
results 

It depends on the characteristic of the medium L and the beam 
input displacement .r, but is independent of wavelength. Let us put 
some typical numbers for a sequence of closely spaced gas lenses. 3 For 

L = \m 

Xi = 2mm 

the range of validity results 

z «4300w. 

The inequality (8) determines the amount of perturbation for which 
the calculations are applicable. For a = 4, that inequality reads 

X 1 
a L*P 



and with the help of (40) 



a\«™. (42) 

L ex { 



The physical significance of ak/L can be derived from (1). With only 
fourth-order aberration the refractive index of the perturbed medium is 



n = n, 



-d~-KT 



The ratio between the perturbing term and the quadratic term at an 
ordinate x is 



'irx 
r x = a\ 
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For x equal to a beam half-width W = VaL/tt, 

» = a|- (44) 

With (42) the following self-explanatory inequality is derived 

r w «— • (45) 

ex { 

Let us continue with the previous example. For 

X = 1m 
and 

W 

r " - "- 1 S • 

we obtain 

r„- = 0.0058. 

The pseudo-period (22) results 



D = |- f (46) 



that is Z) = 460 m. 



VI. CONCLUSIONS 

An off-axis beam oscillates and reproduces itself periodically for any 
arbitrary length only when it propagates in the hyperbolic secant 
square lens-like medium described by A. H. Carter. Nevertheless, 
within some approximations the quadratic law yields similar results 
with much simpler mathematics. For example, a Gaussian beam which 
is properly matched to maintain a constant spot size will oscillate 
periodically with period 2L and amplitude x { . We have found that 
this approximation holds for a distance z that obeys the inequality 

We surmise that the inequality is valid also if the width of the beam 
is not perfectly matched and if the guiding medium is not continuous 
but made of a sequence of square law lenses. For a typical sequence 
of gas lenses with half natural period L = 1 m and beam displacement 
Xi = 2 mm the range of validity is 

z « 4300m. 
A more realistic lens-like medium though, is one in which the quadratic 
law is perturbed by aberration terms of higher order. Again, a Gaussian 
beam of half-width W = VxL/r entering parallel to the axis at a 
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distance x t from it oscillates at the beginning with period 2L, but 
changes shape slowly as it travels along, increasing the width of the 
beam up to 2(x { + W), after which the beamwidth shrinks back to 
the starting value 2W. This process repeats at intervals D that depend 
on the nature of the aberration. For fourth-order aberration 

*-!£■ 

where r w is the ratio between the fourth-order term and the quadratic 
term of the refractive index at a distance from the axis equal to the 
natural half beamwidth IT = va/Vx of the medium. These results 
hold as long as 

A, W 

r w « 

eXi 

Where the beam is large and distorted, the field intensity at one cross- 
section exhibits several maxima. The number of them and their re- 
solution varies along the trajectory. Both increase with the ratio 
of Xi/W. 

We saw that one of the effects of the aberrations is to smear the 
beam size to roughly twice the displacement that the same beam would 
have in an ideal quadratic medium; furthermore, it is known that in 
a sequence of perfect but randomly misaligned lenses the rms deviation 
of a beam from the axis grows proportionally to the square root of 
the number of lenses. 14 Therefore, if the lenses have aberrations we 
conclude tentatively that the rms beam size grows with the same law. 
If redirectors 15 are to be used to compensate for misalignment of the 
lenses the corrections probably must be made before a large break-up 
of the beam occurs. 

APPENDIX 

Approximate Solution of 

t + { P + i-(iy +?fl ..(r-'[/„ -(i)>. -o. 

Provided that 

Y,a a {\/LY ,2 - X f n «p + \ (47) 

and 

E a a (\/L) a/2 -\S/2) a - 2 « 1 (48) 

up to a value of £ to be defined later, the differential equation can 
be solved using the stationary perturbation theory. 12 
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The solution given in terms of parabolic cylinder functions is 

E p = D,(Q+ T,c n D n ®. (49) 



r,Hp 



The summation is a perturbation on D v {£) and it extends from to « 
excluding the pth term. 

The eigenvalue corresponding to the pth eigenfunction E p is 

In order to find c n and /„ we substitute E p given by (49) in the dif- 
ferential equation. Neglecting terms containing products of a a c n and 
knowing that the parabolic cylinder equation is 



at 

we obtain 



d 2 D p (£) + 



> + §-« 



DJ® = 0, (50) 



E (p - n)c n D n ($) + £ a„(|) Q/ '[/„ - (|) B ]^ft) - 0. (51) 

The fimctions D„(£) with integer index q, are orthogonal, therefore, 
multiplying (51) by !><,(£) and integrating from — oo to + °o one derives 



c„ = 2- «-lr ) ^ 

(p - n) / Dl(Q # 

•/ — CD 



(52) 



and 



/ a = ■ / - m y 7 (53) 

/" Z * B * 

After integration 16 the exphcit results of the last two expressions are 
/ x \-/2-. 2"+-— '"qlpl 
? a 'U p-n 



c„ = 



Z 

tn-0 



/ p + n - a ^ (54) 

r- 2 — J ! 



(p — m) ! (w — to) ! to ! ( m 

if p + 7i + a is even, 
.0 if p + to + a is odd 
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and 



U = - 



2~ 3a/2 a\p\ 2^ 7 r~ if a is even, 

'"°(m!) 2 (p -w)!(|- m)\ (55) 

P if a is odd . 



These results are applicable as long as the inequalities (47) and (48) 
are satisfied. Both are harder to satisfy for large values of /„ and £ 
which we proceed to find next. 

According to (55), /„ increases with the order of the mode p and 
for large va ues of p 

U=- J rx 2 v" /2 . (56) 

"■&') 

The function E v given in (49) has significant amplitude essentially 
in the same range of £ than the unperturbed solution D p (£). Since 11 

D&) = (-!)" exp (£) ^ exp (- |) (57) 

a good approximation for £ » 1 is 

Z) p (£)^exp(-£74)r. (58) 

This function has significant amplitude for values of £ smaller than 
that of the second inflexion point. We find it by making the second 
derivative of (58) equal to zero and by finding the largest solution of 
that equation. The result, again for large p is 

Ux = 2Vp. (59) 

Substituting j a and £ maI of (56) and (59) in (47) and (48) we obtain 

Hs«) 

and 

A W2-1 

Za Q (^pJ «1. (61) 
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Furthermore, 

< 1 



r 1 1 



for 



> 0. 



Therefore, condition (61) is the most stringent. Given a a , X, and 
L it establishes which is the highest-order mode p for which the per- 
turbation calculations apply. 
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